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ABSTRACT 

A strong candidate for a source of gravitational waves is a highly magnetised, rapidly 
rotating neutron star (magnetar) deformed by internal magnetic stresses. We calculate 
the mass quadrupole moment by perturbing a zeroth-order hydrostatic equilibrium by 
an axisymmetric magnetic field with a linked poloidal-toroidal structure. In this work, 
we do not require the model star to obey a barotropic equation of state (as a realistic 
neutron star is not barotropic), allowing us to explore the hydromagnctic equilibria 
with fewer constraints. We derive the relation between the ratio of poloidal-to-total 
held energy A and ellipticity e and briefly compare our results to those obtained using 
the barotropic assumption. Then, we present some examples of how our results can 
be applied to astrophysical contexts. First, we show how our formulae, in conjunction 
with current gravitational wave (non-) detections of the Crab pulsar and the Cassiopeia 
A central compact object (Cas A CCO), can be used to constrain the strength of the 
internal toroidal fields of those objects. We find that, for the Crab pulsar (whose 
canonical equatorial dipole field strength, inferred from spin down, is 4 x 10 8 T) to 
emit detectable gravitational radiation, the neutron star must have a strong toroidal 
field component, with maximum internal toroidal field strength B tnl = 7 x 10 12 T; 
for gravitational waves to be detected from the Cas A CCO at 300 Hz, B tnl ~ 10 13 
T, whereas detection at 100 Hz would require B tm ~ 10 14 T. Using our results, we 
also show how the gravitational wave signal emitted by a magnetar immediately after 
its birth (assuming it is born rapidly rotating, with A < 0.2) makes such a newborn 
magnetar a stronger candidate for gravitational wave detection than, for example, an 
SCR giant flare. 

Key words: MHD - stars: magnetic field - stars: interiors - stars: neutron - gravi- 
tational waves 



1 INTRODUCTION 



A str ong magnetic f ield can deform a star (jchandrasekhar fc Fermill953l : lFerrarcll954[GoosenJl972l:lKat jl989l:|Pavne fc Melatosl 



200d; Haskell et al. 2008h . The ellipticity is roughly proportional to the magnetic energy ( Haskell et al. 20081 : 



DalPOsso et al 



200£ 



Therefore, it is expected that magnetars, with their intense magnetic fieldfl possess significant ellipticities. This makes 
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The magn etic field strengths of ma gnet ars are inferred from their periods and spin-down rates, assuming dipole spin-down, but see, 
for example, iKouveliotou et alj il998l ) and iThompson fc Blaesl jl998h for an alternate method, applied to SGR 1806-20. 
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them good candidates for gravitational wave s ources (jBonazzola fc Gourgoulhon 19961 ; Melatos fc Payne] 2005 ; Stella et al 
20051 ; lHaskell et alJlioollDall'Osso et alj|2009h . 



In this paper, we explore the relation between magnetic field configuration (particularly the relative strengths of the 
toroidal and poloidal components) and ellipticity. We do this by starting with a model star in hydrostatic equilibrium and 
impose an axisymmetric magnetic field on it, matched to an external dipolar field to ensure a well-posed problem. We treat 
the magnetic force as a perturbation on this equilibrium and calculate the resulting small changes in pressure and density. 
A similar calculation of the ellipticity has been done for purely toroidal fields in type II superconducting neutron stars by 
Akgiin fc Wasserman (2008). Theoretical arguments, as well as the properties of pulsar glitches, strongly suggest that super- 
fluids (likely including superconducting protons) exist in the interior of neutron stars (Link, Epstein, & Baym 1993; Yakovlev, 
Levenfish, & Shibanov 1999). While superfluidity may have a significant impact on the allowable magnetohydrodynamic equi- 
libria, we ignore this effect in the present study [as also done, e.g., by Haskell et al. (2008), Lander & Jones (2009), and Ciolfi, 
Ferrari, & Gualtieri (2010)]. Also, while our field (described in Sec. 2) superficially resembles the 'twisted torus' of Ciolfi, 
Ferrari, & Gualtieri (2010), our calculation is purely non-relativistic and, unlike Colaiuda et al. (2008), we do not include 
quadrupolar (and higher-order multipoles) fields in our calculation. 

The main novelty of this paper, compared to previous works on the subject [such as Colaiuda et al. (2008), Haskell et 
al. (2008), Lander & Jones (2009), and Ciolfi et al. (2010)], is that we do not require the star to be barotropic, i.e., we do 
not require it to have a one-to-one relation between pressure and density, p — p(p). The reason for this is that neutron star 
matter is a multi-species fluid, in which different kinds of particles co-exist, at the very least neutrons, protons, and electrons. 
The relative abundances of these particles can be adjusted by weak interactions (such as beta decays; Urea processes in 
astrophysical jargon) and diffusive processes. However, the time-scales for these processes are much longer than the time for 
the fluid to reach a hydromagnetic equilibrium (force balance) state, which is likely to happen in a few Alfven times (of order 
seconds). In the process of settling to this hydromagnetic equilibrium, each fluid element is free to move following the relevant 
forces, but it conserves its composition (relative abundances of species) , which will thus not be exactly a function of pressure 
and density. Once the equilibrium is reached, the star will be radially stratified 1 Pethick 19921 ; Reisenegger fc Goldreich 1992 *1 
to zeroth order, but the spherical symmetry (and with it the stratification) will be slightly perturbed by the presence of the 
magnetic field. 

On time-scales intermediate between the Alfven time and the weak interaction and diffusion time-scales, the star will thus 
be in a hydromagnetic equilibrium state in which the composition is not strictly determined by the density or pressure, and 
the latter thus do not have a one-to-one correspondence with each other. The present study applies to this regime, considering 
a hydromagnetic equilibrium state described as a spherically symmetric, non-magnetic background in which both pressure and 
density depend only on the radial coordinate, plus small, non-axisymmetric pressure and density perturbations, caused by the 
presence required to balance the Lorentz force, and which are not directly related to each other, due to the position-dependent 
composition. 

In this state, the particles will not be in equilibrium with respect to weak interactions or diffusive processes. The latter 
will act on longer time-scales, changing the c omposition and thus the density and pressure per turbations, and thus causing 
a long-term evolution of the magnetic field ( Goldreich fc Reisenegger] 19921; iReiseneggerl I2OO9I) . The s e pro cesses have been 
studied in some detail in a simplified, one-dimensional model ( Hovos. Reisenegger. fc Valdivial [20081 . 2010l ). which is being 



extended to a more realistic, axially symmetric geometry (jBecker et al.ll201lh . When these non-ideal MHD processes become 



relevant, the ideal, non-barotropic MHD equilibrium studied in this paper most likely breaks down. At least during an 
im portant part of this later evolution, the ne utron star matter is likely to be describab l e as t wo independent fluids, as done 
bv lGlampedakis. Andersson. fc Landed ( 2011 ) and lLander. Andersson. fc Glampedakisl 1 2011 ). 

Analytically, imposing barotropy means requiring that th e magnetic force per unit mass is expressible as a gradient 
I PrendergastHl95a : Monaghanll 19651 : Akgiin fc Wass ermarfeOQS ) . Rela xing this constraint for the reasons explained above 
enables us to explore a greater range of equilibria ([Reisenegger 20091 ) Fl 

In Section 2, we describe our model in detail, starting with the zeroth-order hydrostatic equilibrium and perturbing it 
magnetically. In Section 3, we derive a simple formula relating the ellipticity e to the magnetic energy, the poloidal-to-total 
magnetic energy ratio, and the mass and radius of the star, for steady states with a par abolic density profile or a polytropic 



equation of state with index n = 1. We also compare our results to those obtained by lHaskell et al] (|2008t ). who assumed 
barotropy. The formula for e has several astr ophysical appl ications. In Section 4, we discuss these and compare our work 
to previous studies in the literature, such as ICutlerl (|2002h . We apply our results to gravitational wave emission from a 
newly born, rapidly rotating magnetar. We discuss the possibility of detecting the gravitational-wave signal from a newly 
born magnetar in the Virgo cluster and relate its detectability (quantified by the signal-to-noise ratio) to the magnetic field 
configuration, emphasizing the poloidal-toroidal energy ratio. We then apply our formulae, together with the current upper 



2 Ther e are other mechanisms which can defo r m a neutron star, such as 'mountain s' on the elastic crust, formed by wav y electron-capture 
layer s l| Ushomi rskv. Cutler, fc BildsteJ l200d: lHaskell. Jones, fc Anderssorj I2OO6I ) . rotation, or crustal shear stresses l|Pines fc Shahaml 
H;[Cn 



In this paper, we ignore these other mechanisms, so as to concentrate on the magnetic deformation. 



Gravitational radiation from a magnetically deformed non-barotropic neutron star 3 



limits from the Laser Interferometer Gravitational- Wave Observatory (LIGO), to constrain the strength of the toroidal fields 
of the Crab pulsar and the Cassiopeia A central compact object. We summarize our results and discuss further applications 
and refinements of our idealised model in Section 5. 



2 HYDROMAGNETIC PERTURBATION 



Let ( r, 9, <j>) be the spheri cal polar coordinates, with r expressed in units of the stellar radius R, (so that it is dimensionless) . 
As in lAkgun et al.l (|201l|), we choose a magnetic field configuration such that: (1) V • B = is satisfied; (2) the magnetic field 
is axially symmetric around the z-axis; (3) the poloidal component of the field is continuous with an external dipole field, 
which vanishes as r — ► oo; (4) there are no surface currents (the tangential, toroidal component vanishes at the surface); (5) 
the toroidal component is confined to the region of closed poloidal field lines around the neutral line (on which the poloidal 
component vanishes); and (6) the current den sity remains fini t e and continuous everywhere in the star. We write the magnetic 
field inside the star in the form pioneered by Chandrasekhar ( 19561 ). as a sum of toroidal and poloidal components, 



B = Bo[r]pVa(r, 9) x \7cj> + rjt/3(a.)\/4>], (1) 

where r\ v and r) t are dimensionless parameters which define the relative strengths of the poloidal and toroidal components, 
respectively. The flux function a(r,9) is taken to be f(r) sin 2 9. Note that this particular form of a(r,9) is only applicable 
when we try to match our field to an external dipole field; other multipoles will have different ^-dependences. The radial 
dependence of a is defined by 



fir) 



35 



6r 4 3r' 



(2) 



The function f(r) is postulated to be of this form to ensure that the field described by Eqs. (l)-(3) is continuous with a 
dipole field outside the s tar, that ther e are n o surface currents, and that the current density is finite at the origin [for a more 
thorough derivation, see lAkgiin et all |201lh ]. The function j3(a) is chosen as 



(a 



0(a) = ^ 







a < 1, 



(3) 



which ensures that the toroidal field is confined to the region where a exceeds 1, the value taken by a at r — 1, 9 — it/2, with 
the current density going continuously to zero at this boundary. 

The magnetic structure described by Eqs. (l)-(3) is drawn schematically in Fig.[T] The field lines lie on surfaces of constant 
a. The toroidal component is completely restricted to the shaded region, which is bounded by the surface a(r,9) = 1. Note 
that both the toroidal and poloidal components are continuous across a = 1, so that there are no surface currents at a — 1, 
either inside the star or on the surface at the equator. The poloidal component inside the a 1 region (the shaded region) 
lies on nested tori around the neutral line, with the strength of the poloidal component continuously decreasing towards the 
neutral line. Outside the star, the poloidal component is matched to a dipole field, consistent with a(r, 9) oc sin 2 9. 

If the magnetic field is treated as a perturbation on the background^, we can write the force balance equation to first 
order in B 2 / pop as[3 



(l//xo)(V x B) x B = VSp + 6p\7<S>, (4) 

where the Cowling approximation has been taken (i.e., 5& — 0). Note that we do not require the density perturbation Sp to 
be a function of the pressure perturbation Sp (the barotropic assumption). Therefore, the magnetic field is not required to 
satisfy the Grad-Shafranov equation. 

Substituting the field given by Eqs. (l)-(3) into Eq. (4) gives 



3 Care must be taken, however, that neither the poloidal nor the toroidal component of the field approaches the virial limit, 

where the magnetic field energy is equal to the gravitationa l binding energy jLattimer & Prakash 2007). The exact value of this limit 
depends on actual field structure. In fact, [Rciscncggcr ( 2009|) showed that the upper limit is set by stable stratification at ~ 10 13 T. 

4 Throughout this paper, we use SI units. 1 T = 10 4 G. 
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por 2 sin 2 



-(r/pVaAa + % 2 /3V/3) = VSp + SpV$, (5) 



with the operator 



Q r 2 r 2 QQ ^0 QQ ^ 

The ^-component of the right-hand side of the force balance equation (5) gives d(Sp)/d8, which is integrated with respect to 
8 to give Sp. Then Sp is obtained from the radial component of (5). Thus, we write 



,__9 2 _ sing d f 1 d . 



dp— = 



-Bja 
por 2 sin 6 



77 2 Aa + 2r/ 2 ( y - ^- + 3a - In a + C 



(7) 



(8) 



where Co is a constant, needed to ensure Sp is continuous across a = 1. In this particular case, where the field is defined by 
Eqs. (l)-(3), we find Co = —11/6. Equations (7)-(8) are applicable to the region where a ^ 1; to obtain similar expressions 
for the a < 1 region (where the toroidal field vanishes), one simply sets r)t = 0. 

An arbitrary function of radius h(r) can be added to 8p. This is equivalent to changing the background profile by a 
spherically symmetric function, hence it does not affect the mass quadrupole moment, which is the focus of this paper. 
For computational simplicity, this function is set to zero (so that Sp vanishes along the z-axis). Note that, while both the 
Lagrangian and Eulerian pressure perturbations vanish at the surface, the Eulerian density perturbation [Eq. (8)] does not 
(only the Lagrangian perturbations need to be zero at the perturbed surface). 

We integrate the squares of the toroidal and poloidal components of the magnetic field over all space (including the 
vacuum, to r -> oo) to obtain their energies. We can write the total magnetic energy E m as 



, , 2 , 2-,^-Bo-R* 

E m = (<V7 P + a t r) t ) , (9) 

po 

with the dimensionless constants a p = 8.48 and a t = 1.65 x 10~ 5 . We can then express the ratio of the energy of the poloidal 
magnetic field component to the total magnetic energy as A 

A= J- y «*rB;/*, = _gg_ 

J v dVByp vi+qvV 1 ' 

with q = 1.95 x 10 -6 . By definition, one has ^ A ^ 1, where A = and A = 1 represent a star with a purely toroidal and a 
purely poloidal field, respectively. 

In Fig. [21 we plot Sp and 5p(d&/dr) for various values of A, in units of B 2 /po, taking rj p = 1 without loss of generality. 
In the first column, we plot pqSp/Bq (solid curves) and po5p(d&/dr)/ B$ (dashed curves) along 8 = ty/2 as a function of r for 
some values of A, to show how the internal toroidal magnetic field gradually changes the pressure and density distributions. In 
the purely poloidal case, Sp has a maximum at r = 0.54. As the toroidal field strength increases (i.e., as A < 1), local minima 
for Sp and Sp(d&/dr) develop straight away at r — 0.78 (where the toroidal part of the magnetic field has a maximum), which 
eventually become global minima. In a sense, the pressure and density perturbations which are solely due to the poloidal field 
(the peak at r — 0.54) are gradually overwhelmed by the perturbations caused by the strengthening toroidal field component 
(the peak at r = 0.78). In the second and third columns, we plot contours of poSp/Bg and po5p(d&/dr)/ B$ in a meridional 
quadrant, to show their dependence on the polar angle. Note how the contours are clustered around the new minimum at 
r = 0.78, created by the increasingly stronger toroidal component, as A decreases. 



3 MAGNETICALLY INDUCED DEFORMATION 
3.1 Ellipticity 

The gravitational wave strain generated by a rotating rigid body is directly proportional to its ellipticity. Hence, our first step 
is to calculate the ellipticity of the star due to the magnetic configuration given by Eqs. (l)-(3). The ellipticity e is defined as 
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Figure 1. A diagram of the magnetic field lines described by Eqs. (4)-(6). The surface of the star is represented by the dashed circle. 
The toroidal magnetic component is confined to the shaded region and describes a circle around the z-axis. The poloidal field lines vanish 
at the neutral line. All field lines lie on surfaces of constant a. 



where 7o is the moment of inertia of the spherical star. The moment-of-inertia tensor is given by 

I jk = R* [ dVlp^ + Sp^e^djk-xjXk), (12) 
Jv 

which allows us to rewrite (11) as 

e = irR 5 Jd 1 / drd6Sp(r,6)r 4 sin6(l - 3cos 2 6>). (13) 
Jv 

Roughly speaking, the ellipticity induced by the magnetic field is proportional to the magnetic energy, which determines 
<5p, as is evident from Eq. (8). We wish to know how the details of the magnetic configuration, specifically A, control e. 
Below, we apply our model to two different unmagnetized steady states. We show the density and d&/dr profiles for these 
unmagnetized steady states in Fig. [3] (a) and (b), respectively. 



3.2 Parabolic steady state 

In the absence of a magnetic field, the star is in a spherically symmetric hydrostatic equilibrium. Here, we consider the simple 
density profile l|Akgun et alJl201ll ) 

p = p c (l-r 2 ), (14) 

where p c = 15M*/(8ttRI) is the density at the centre. The gravitational acceleration (i.e., the gradient of the gravitational 
potential V"! 1 ) corresponding to this steady state is 

where M» and R, denote the mass and radius of the star, respectively and G is Newton's gravitational constant. From the 
gravitational acceleration, we obtain the pressure profile 
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Figure 2. (Left column) Pressure and density perturbations fiQ&p/ Bq (solid curve) and [io8p\d<& /dr\ / Bq (dashed curve) versus normalised 
radius r, at 9 = n/2, with r) p = 1 and various values of A. (Middle and right columns) Contour plots of [iq&p/Bq and fioSp\d&/dr\/ Bq, 
respectively, on the meridional quadrant. Note that contour levels differ between panels, they are auto-scaled to emphasize the structure 
of /iqSp/Bq and poS p\d& / dr\ / Bq . The dashed curves represent the surface of the star. 
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Figure 3. (a) Density profiles for the unmagnctized parabolic steady state given by Eq. (14) (solid curve) and the unmagnetized n = 1 
polytropic steady state given by Eq. (18) (dashed curve) in units of M t /R^; (b) gravitational acceleration d^/dr as a function of r for 
the unmagnetized parabolic steady state given by Eq. (15) (solid curve) and the unmagnetized n = 1 polytropic steady state given by 
Eq. (20) (dashed curve) in units ofGAh/Rl. 



P = Pc(l-^ + 2r 4 -0, (16) 

where p c = /{16nRi) is the pressure at the centre. We emphasize that this is a particular, simple choice of density 

profile, chosen to render the following calculations tractable; it is not motivated directly by physical arguments or observations. 
For this steady state, Eq. (13) gives 

e = 5.98 x ID" 6 ( V (_^_) " ( M 4 (17) 

V5xl0 10 T7 \1AM Q J V 1 o 4 my V A / 

where -Bsurface is the surface magnetic field strength at the equator; we have -Bsurface = rj p Bo in our formulation. 



3.3 Polytropic steady state 

For comparison, we also consider the n = 1 polytropi c star, whose zeroth-order density and pressure configurations are given 



by l|Akgun fc Wassermanll2008l ; iDall'Qsso et al 



ytropic 
. 2009) 



p c sin(Trr) nM, 
P= nr ' Pc = W (18) 

p = kp 2 , = (19) 



The gradient of the gravitational potential is 

d$ GM, 

For this steady state, Eq. (13) gives 



dr 7rJ?* 



(20) 



Bsurfece \ { M, \ ( R* Yf 1 0.385 \ ( . )[} 



5xl0 10 T/ \lAM e J V 104m / V A 
We plot e as a function of A for an n = 1 polytropic star as the dashed curve in Fig. [4] Note that, like in the parabolic 
steady state, the ellipticity in the n = 1 polytropic steady state case also vanishes at A c ~ 0.39. Equations (17) and (21) 
show that the n = 1 polytrope and the parabolic steady states give ellipticities within 5% of each other. For consistency and 
analytic simplicity, throughout the rest of this paper, we use the parabolic steady state to calculate ellipticity. 
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Figure 4. Ellipticity e (in units of 10 6 ) versus A, the ratio of poloidal magnetic field energy to total magnetic field energy, for r] p = 1, 



B = 5 X 10 10 T, R, = 10 4 m, M* = 1.4M . Note that e is proportional to BgR 4 M»~ 
[Eq. (15)] and n = 1 polytrope [Eq. (19)] steady-state density profile, respectively. 



The solid and dashed curves are for a parabolic 



3.4 Comparison with Haske ll et alj (2008) 



We now briefly compare our result to that obtained by Haskell et al. 1 20081 ). The latter authors postulate an 



1 polytropic 

star which also obeys the barotropic condition. This restricts the form of the magnetic field, as the (f i rst-or der) density and 
pressure perturbations are required to be proportional to each other. The magnetic field in Haskell et al. (2008) is characterized 
by the discrete eige nvalue A (not to be confused with our A), as described by Eqs. (B.6)-(B.10) of that paper. Furthermore, 
Haskell et al.l (|2008r ) do not make the Cowling approximation. Otherwise, the two methods are fundamentally similar: begin 
with an unmagnetized steady state in hydrostatic equilibrium, apply an axisymmetric magnetic field, assume the departure 
from the hydrostatic equilibrium caused by the magnetic field is small enough to be treatable as a perturbation, then calculate 
the perturbations to pre ssure and density due to the applied magnetic field. 



Haskell et al 



1 200&] ) present their results for e as a function of A (and, hence, volume-averag ed absolute values of the 
toroidal and poloidal field strengths B t and B p ) in their Table 1 [see also the associated Erratum (lHaskell et alJboogl )]. We 
calculate e for our field configuration using Eq. (21) and compare our e to those obtained by lHaskell et al] |2008l ) for the 
same values of internal poloidal-to-total field energy ratios in Tab le 1 (where we also show the corresponding values of A). 
We emphasize here that a direct comparison between the results of Haskell et al. ( 20081 ) and ours is not entirely appropriate, 
because we use very different field configurations. However, by using Table 1, we simply wish to illustrate how the resulting 
elliptici ties are very d i fferen t, even when using the same poloidal-to-total field energy ratios. 

As Haskell et al. 1 20081 ) did, we find that the star is always prolate for the chosen field configuration and toroidal field 
strengths, and becomes more prolate as the toroidal field component gro ws. We general l y find larger values of e for smaller 
A and smaller e for larger A. We caution the reader again that we and lHaskell et al. (2008) use different magnetic field 
configurations: 



• our toroidal field is concentrated within an equatorial torus, while that of Haskell et al. 1 20081 ) is not confined to this 
region, varies smoothly across the stellar interior (with the number of oscillations controlled by A), and vanishes at the origin, 
the stellar surface, and the z-axis; 

• our poloidal field is continuous with an external dipole field, whereas that of lHaskell et al. I |200Sl ) vanishes at the surface; 

• the comparisons we make in this section are between states with the same internal poloidal-to-total field energy ratios. 



While our methods also differ in that we take the Cowling approximation whereas Haskell et al. 1 20081 ) do not, we contend 
that it is the fundamental difference between field configurations that causes the discrepancies in e seen in Table 1. To 
illustrate this further, we show the contour plot of (j,o5p(d<&/dr) [i.e., the radial component of the force balance equation Eq. 
(4)1 for A = 7 . 459 a s Fig. [5] Note how the radial component of the internal poloidal field lines of the configuration used by 
Haskell et al. (20085) vanishes at the surface of the star (Fig. 0c, left panel), unlike our field configuration, where the radial 
component is continuous with a dipole field outside the star (Fig. [5]c, right panel). As an aside, in the context of accreting 
neutron stars, Haskell et al. (2006) have shown that abandoning the Cowling approximation alters the possible maximum 
surface deformation (due to accreted mountains), and ellipticity, by a factor of 0.5-3 (depending on the density profile chosen). 
While this is a different physical phenomenon than the one discussed here (i.e., deformation due to accretion rather than 
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Table 1. Comparison between the ellipticities e obtained by Has kell et all l|2008h (eh; fourth column) and by using Eq. (21) (e2i, 
fifth column), as a function of the eigenvalue A and volume-averaged toroidal field strength Bt, with the volume-averaged poloidal 
field strength B p kept constant at 10 8 T. We show the corresponding poloidal-to-total field energy ratio A in the third column to 
make contact with the rest of the calculations and results in this paper. 
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-2 x 10" 
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due to internal magnetic fields), this result indicates that taking the Cowling approximation may change our ellipticities by 
a factor of order unity. Of course, a full calculation needs to be undertaken to determine if this is truly the case, and, if it is 
indeed so, whether the Cowling approximation increases or lowers ellipticities. 



4 GRAVITATIONAL- WAVE ASTROPHYSICAL LIMITS ON e AND A 

Equations (10) and (17) allow us to use observational limits on ellipticity to constrain the strength of the internal toroidal 
field component. There are many possible contexts in which we can apply these ideas and results. In this section, we look at 
two illustrative applications: 

• isolated pulsars with or without timing solutions; 

• birth of isolated magnetars. 



4.1 Upper limit of gravitational radiation from the Crab pulsar 

Recent Laser Interfero meter Gravitational- Wave Observatory (LIGO) Science Run 5 observations of the Crab pulsar have put 
limits on its ellipticity (jAbbott et aljboostliuioh . Using Eq. (17) and these data, we can constrain the strength of the surface 
and internal magnetic fields of the Crab pulsar. In Table 2, we list the lower limits of A based on t he late st upper limits o f 
detection from LIGO and the planned Einstein Telescope in its B ( Hild. Chelkowski. &: Freisei 20081 ) or C (|Hild et al-lboid ) 
configuration. From the fundamental definition of A [Eq. (10)] and the magnetic field defined by Eqs. (l)-(3), we can express 
the maximum of the internal toroidal field B tm as a function of Bq and A to facilitate comparison: 



B tm = 18.2B [t-A ■ ( 22 ) 

To give a numerical example, for B S urfaco (the surface dipole field at the equator, equivalent to r/pBo in our formulation) 
corresponding to the theoretical value for pure dipole braking (namely 3.8 x 10 8 T) and A = 0.4 (i.e., Btm = 8.4 x 10 9 T) 
we find e = 8.7 x 10~ 12 , which is much lower than the LIGO upper limit. It is clear that, if indeed the surface magnetic field 
strength is -B sur f acc ~ 10 8 T, as inferred from electromagnetic spin down, detection of gravitational waves from the Crab pulsar 
is only possible if the star is prolate with a very strong toroidal component (A < 10~ 6 , corresponding to Btm > 1-82 x 10 12 
T). In other words, the poloidal field alone cannot deform the Crab pulsar enough to be detectable as a gravitational wave 
source by LIGO and the planned Einstein B or C; a strong internal toroidal field is required to deform the star into a high-|e| 
prolate shape. Indeed, if gravitational waves were to be detected from the Crab pulsar, the analysis presented in this paper 
offers a method to infer the strength of this internal toroidal field. 



4.2 The case of Cassiopeia A Central Compact Object 

The supernova t hat gave birth to the supernova remnant Cassiopeia A (henceforth Cas A) was probably observed in 1680 
( Ashworthlll980h . although the central compact object (CCO), hy pothesized to be a neutron star, was only positively iden- 
tified as a bright X-ray po int source in 1999 [jTananbaum| [l999). making this one of the youngest neutron stars in our 
Galaxy (|Ho fc Heinkelboogl ). Because the location of the CCO is well known, it makes for an attractive target for a directed 
gravitational wave search with LIGO ( Wette et al. 200ct ). 




Figure 5. Comparison of the field configuration of lHaskell et al,l (|200Sf ) and that of this paper, taking A = 1.77 X 10 2 [corresponding 
to A = 7.459 in the notation of Hask ell et al. (2008)]: (a) fi S p(d<S> / 'dr) / 'b\ as a function of radiu s r (normalised to th e stellar radius), at 
the equator for the fields described by Eqs. (4.1)-(4.3) (solid curve) and by Eqs. (B6)-(B10) of Haske ll et al] l|2008h : (b) contour plots 
of p,o5p(d$/dr)/ B t for the field used by lHaskell et all J200ST) (left panel) and our field (right panel); (c) internal poloidal field lines for 
the field used bv lHaskell et al.1 120081 (left naneH and our field fright nanelV 
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Table 2. Inferred minimum A m i n of the Crab pulsar, using Eg, (17), based on th e upper limits of gravitational wave strain (ho) 
from various gravitational wave detectors jAbbott et al.ll2008l , l20icl; lBennetdl2010T) . We note that, in the case of the Crab pulsar, 
because of its (relatively) weak magnetic field strength, significant wave strain can only be achieved if the star is highly prolate. 
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While the location of the CCO is well known, it has no detectable pulse in any electro magnetic band, r endering it 
impossible to measure its rotation frequencjH, frequency derivative, and magnetic field a priori ( Wette et al. 2008). However, 
recen t analysis has shown that a neutron star model with a low magnetic field can explain the Cas A observations (|Ho fc Heinke 
2009). After fitting several different atmospheric models (namely blackbody, pure hydrogen, helium, carbon, and nitrogen) 



to the measured X-ray spectrum, Ho and Heinke (2009) concluded that the Cas A CCO is most likely a neutron star with a 
carbon atmosphere, radius 8-17 km, mass 1.5-2.4 M©, and upper limit of 10 7 T for the magnetic field. For this set of physical 
parameters, our model predicts e ~ 10~ 15 if the internal field is purely poloidal. Therefore, we expect th e star to be p rolate 
wit h a large toro idal component (a very low A) if it is to approach the LIGO upper limits estimated by Wette et al. ( 20081) 
and lWettei(|2010h . 

Based on the Ho and Heinke (2009) conclusions and the most recent uppe r limits on the ellipticity set by LIGO, obtained 
from indirect upper limits on intrinsic gravitational wave strain (|Wettell201oT ). we can now set lower limits on A of the Cas 
A CCO (i.e., upper limits on the toroidal field strength). Assuming M* = 1.5-Mq, R = 10 4 m, we take the upper limits on 
e for four different putative gravitational wave frequencies / (twice the spin frequency): 100, 150, 200, and 300 Hz. We then 
plot the lower limit on A as a function of the (unknown) surface field strength in Fig. [6] The curves in Fig. [6] thus indicate 
the A needed at various B sur f acc for detection by LIGO at 100 (dotted curve), 150 (dashed-dotted curve), 200 (dashed curve), 
and 300 Hz (solid curve); the region to the left and above each curve generates e that is below the s ensitivi t y lim its of LIGO. 
The LIGO upper limits for e at lower / are higher [e.g., e = 3.6 x 10~ 4 at / = 100 Hz, see Fig. 4 of lWettel (feoicl )]. hence the 
smaller limits on A at lower frequencies seen in Fig. [6] 

While the corresponding limits for the internal magnetic field strength may seem to be very large (even approaching 
Byiriai at low /), note that, given the lack of other information to constrain Btm, these are merely upper limits. From Fig. [6] 
and Eq. (22), we see that our model predicts that, for gravitational waves to be detected from Cas A at 300 Hz, the CCO 
must have B tm > 4 x 10 13 T (if it has mass 1.5 M Q and radius 10 4 m) or _B tm ^ 2 x 10 13 T (if it has mass 2.4 M Q and radius 
1.7 x 10 4 m); for gravitational wave detection at 100 Hz, the CCO must have Btm 1.2 x 10 14 T (if it has mass 1.5 Mq and 
radius 10 4 m) or B tm 5? 6.5 x 10 13 T (if it has mass 2.4 M Q and radius 1.7 x 10 4 m). 



4.3 Detectability of gravitational wave signals from newly-formed magnetars in the Virgo cluster 



Acco r ding to certain evolutionary scenarios, magne tars are hypothesized to be born with periods P ~ 1 ms ( Thompson fc Duncan! 
19931 : lYi fc Blackmanl ll998: Dall'Os so et al.ll2009f). If so, significant amounts of gravitational radiation may be emitted dur- 



ing a magnetar's early days ( DaH'Osso fc Stella 2007 : Dall'Osso et al. 20091 ). Gravitational radiation is maximal when the 
symmetry axis is orthogonal to the spin a xis, and the spin down of a neutron star (including a magnetic contribution with 
arbitrary braking index n) is then given by (|Shapiro fc Teukolskv|[l983l : lMelatoslll997l : lwette et al.ll2008l : lDairOsso et al.ll2009l : 
Chung et al.ll201ol ) 



r\n 3 r»6 r?2 
p _ z ™ -"surface 

3/x Ic 3 P n ~ 2 

and the rate of rotational energy loss is given by 
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where, again, Surface is the s urface dipole field at th e equator (equivalent to i] p Bq), and n is the braking constant. 



(24) 



It has been suggested bv lThompson et al.l (2004) that a newly born magnetar can be a good candidate for gravitational 
wave detection during the first 100 s of its life. Firstly, using our model, assuming that the magnetar is born with initial 
period of 1 ms, surface equatorial dipole field of 5 x 10 10 T, radius of 10 km, mass of 1.4 M Q , A = 0.80 we find that the effects 



5 However, due to the matched filtering method used by LIGO, detection of gravitation al waves would a utomatically give rotation rate. 

6 This value is the upper limit of A for a stable star with axisymmetric magnetic field (Braithwaitc 200^). 
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^surface/10 T 

Figure 6. Contours of constant gravitational wave frequencies plotted on the A— -B sur f ace plane, indicating the upper limits on B S urface 
and lower limits on A for each given frequency (solid curve: 300 Hz; dashed cu rve: 200 Hz; dashed-dotted curve: 150 Hz; dotted curve: 
100 Hz) for the Cas A CCO, inferred from the most recent LIGO upper limits llWettell2O10l l and the analysis by Ho and Heinke (2009). 
We take M* = 1.5M©, R = 10 4 m. For the case of M* = 2.4M Q , R = 1.7 X 10 4 m, A is higher by a factor of 3.26 for each given B p and 
gravitational wave frequency. 



of gravitational spin down on the period are negligible compared to the electromagnetic spin down. Spin down purely due to 
gravitational wave emission can only increase the newborn magnetar's period by 10 -3 %, wh ereas pure electromagne tic spin 



down can more than double the period within these first 100 s [cf. top left panel of Fig. 5 of iThompson et"aL ( 20041 )] . Even 



though the effect on spin down is marginal, the power emitted in gravitational radiation during these first 100 s may still be 
significant. Thus, we plot the gravitational rotational e nergy extracted f r om th is newly born, rapidly rotating magnetar in 
the first 100 s as Fig.[7][cf. bottom left panel of Fig. 5 of lThompson et al.l (|2004T )]. Note the close quantitative correspondence 
between our strongly poloidal case (0.9 < A < 1) and their Case C. Our Fig.Qalso shows that energy output cannot distinguish 
between the 0.9 ^ A ^ 1 case and 0.22 ^ A ^ 0.25. We also plot th e gravitational wave energy upper limits from SGR 1806-20 
giant flare observations and SGR 1900+14 simulated storm search ( Kalmus et alJ l2009) for rule-of-thumb comparison (despite 
their astrophysically different natures). It is clear that a newly-born, fast-rotating magnetar with A 0.2 can be a significant 
source of gravitational waves, stronger than an SGR giant flare, and thus a promising target for gravitational wave searches 

duri ng its f i rst 10 s. 

Cutlerl (|2002l ) and Istella et all |2005l ) [further refined by iDall'Osso etafl (|2009r i] approximated the optimal matchec - 



filter signal-to-noise (S/N) ratio of newly born magnetars in the Virgo cluster by integrating the instantaneous signal strain 
amplitude over the frequency range ft ^ f ^ ff {ft is the initial, 'birth' spin frequency of the magnetar and ff is the final 
frequency, when we assume the most significant gravitational wave emission ceases) and averaging over detector antenna 
pattern. For Advanced LIGO at frequencies 0.5-2 kHz, they give the following formula for S/N: 



S/N=l 



ttGI 



_ (Kgw\ 

6c 3 DS X,2 \ Kd ) 



1/2 



21nA 
If 



In 



a + ff 



1/2 



(25) 



where So ~ 2.1 x 10 53 Hz 1 is the squared noise spectral density at frequencies between 0.5 and 2 kHz ( DaH'Osso et al. 20091 ) . 
Kd = 2-7T 3 ^?»Bg Urface /(3/xo/c 3 ), Kgw = 327r 4 Ge 2 //(5c 5 ), a = Kd/ir 2 Kqw , and D is the distance to the source. 

Taking D = 20 Mpc, ft = 1.03 x 10 3 Hz, f f = 0.1 Hz, M« = 1.4 M , R, = 10 4 m as an example, also assuming the 
source is an orthogonal rotator, we plot A as a function of Bsurface [using Eq. (17) to relate A a nd e] for S/N = .3, 1, 3, and 
10 in Fig. El Note that we obtain generally smalle r S/N compared to the predictions made by DaH'Osso et al. 1 20091 ). This 
is because Dall'Osso et al. (2009) used the ICutlerl (|2002l ) formula, which assumes that the toroidal field dominates and the 
poloidal field is negligible, to calculate their e, whereas we use Eq. (17) to calculate ours, which is derived self-consistently 
in Sections 2 and 3 (without neglecting the poloidal component, which makes our star less prolate). According to our model, 
significant detectability (S/N w 10) for Surface = 5x 10 10 T is only possible when A < 0.01, when the star is prolate and the 
maximum toroidal field strength is ~ 3500 times the poloidal field strength at the surface; when the field is purely poloidal 
(A = 1), we predict S/N » 0.18 only. 



Gravitational radiation from a magnetically deformed non-barotropic neutron star 13 



1000 



* 




0.001 



1.00 



Figure 7. Gravitational wave spin-down energy extracted from a rapidly rotating, newly born magnetar (initial period 1 ms, surface 
dipole field 5 X 10 10 T at the equator, radius 10 km, mass 1.4 Mq, distance 10 Mpc) in the first 100 s as a function of the poloidal- 
to-total-magnetic energy ratio A. Also shown for comparison are the gravitational wave energy upper limits from the SGR 1806-20 
giant flare (dotted line: white noise burst waveform 100 ms 100-1000 Hz; dashed line: circularly polarised waveform 200 ms 2590 Hz) 
llAbbott et al. 2008) and from the simulated SGR 1900+14 storm search (dashed-dotted line: circularly polarised waveform 200 ms 1590 
Hz) jKalmus et alj|200sl) . The latter sources are astrophysically unrelated to a newly born magnetar; they are included as rough guides 



to the sort of signal lev els detectable at present. Note that we show the vertical axis in units of 10 51 ergs to facilitate comparison with 
iThompson et al.N2004h . 
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Figure 8. Relation between A and -B sur face for given values of signal-to-noise ratio (S/N), for a M* = 1.4 Mq, R, = 10 4 m source 
located 20 Mpc away, rotating with initial spin period 0.97 ms and final spin period 10 s. Solid curve: S/N = 10; dashed curve: S/N = 3; 
dashed-dotted curve: S/N = 1; dotted curve: S/N = 0.3. 



5 CONCLUSIONS 



Neutron stars in general are not expected to be barotropic 1 Pethick 1992 ; Reisenegger fc Goldreichl 1992). They posses s 



observable external dipole fields and internal linked poloidal-toroidal fields ( Braithwaite fc Nordlundl 



2006; Braithwait 



2009). 



Non-barotropic assumption allows us freedom in constructing such equilibriaQ In this paper, we model a highly magnetised 
neutron star as a fluid sphere in hydrostatic equilibrium perturbed by an axisymmetric, poloidal-toroidal magnetic field, 
whose toroidal component is entirely confined within the star and whose poloidal component i s continuous with a dipole field 



outside of the star. Our model differs from those of previous authors [e.g., Haskell et al. ( 20081 )] in that we do not assume the 



star to be barotropic. We then self-consistently calculate the ellipticity caused by density perturbations in a magnetar and 



7 The stability of these equilibria is not considered in this paper; it is the subject of another paper in preparation jAkgiin et alj|201lf) . 
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derive a simple relation between ellipticity e and the ratio of poloidal to total magnetic field energy A. In principle, equation 
(17) allows one to infer the internal magnetic structure of a neutron star from observable quantities (namely, ellipticity, mass, 
surface field strength, and radius). For a magnetar with surface dipole field strength of -B sur f acc = 5 x 10 10 T, radius 10 km, 
mass 1.4 Mq, with a purely poloidal field configuration (A = 1), we predict an ellipticity of 3 x 10~ 6 and expect such a star 
to be oblate. 

Applying our results to specific situations, we show how current gravitational-wave non-detections can be used to put 
constraints on the strengths of the internal toroidal magnetic fields of the Crab pulsar and the Cas A CCO. Then, we show 
that a newly born, fast-rotating magnetar can emit gravitational radiations with energies of ~ 10 45 J (if A = 1 or 0.2) during 
the first 100 s of its life (on a par with the SGR 1806-20 giant flare) and that a newborn millisecond magnetar in the Virgo 
cluster can be a stronger candidate for a detectable gravitational wave source than an SGR giant flare. We also predict the 
signal-to-noise ratio for a source located 20 Mpc away and relate this ratio to A. 

In this paper, we ignore physical effects such as crustal stiffness, the presence of superfluids, and rotational deformations. 
These effects warrant further investigation. Furthermore , a ful le r model is required to explain d etails of magnetic reconfig- 
uration. In conjunction with works such as Braithwaite ( 20091 ) . Braithwaite fe Nordlund (2006), and iBraithwaite fc Spruitl 
I 2004l . l2006h . our work can also be extended to investigate the implications of magnetic field configurations on various internal 
aspects of a neutron star, such as the equation of state. 
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